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حکیده 


تعیین توزیع مکانی تبخیر و تعرق مرجع (ETO)‏ در مدیریت و برنامه‌ریزی منابع آب» مطالعات 
بیلان آبی و برآورد نیاز آبی گیاهان از اهمیت زیادی برخوردارست. بنابراین هدف پژوهش حاضر 
تعیین مناسب‌ترین روش درونیابی برای برآورد مکانی ETO‏ و پهنه‌بندی این متغیر در استان 
خوزستان است. بدین منظور از داده‌های هواشناسی ۴۲ ایستگاه در یک دوره آماری YA‏ ساله استفاده 
شد و برای تخمین مقادیر ETO‏ در ایستگاه‌ها؛ روش فائو پنمن مانتیث بکار برده شد. به منظور 
برآورد مکانی ETO‏ هفت روش درونیابی مشتمل بر وزنی عکس فاصله. اسپلاین» گرادیان خطی 
سه بعدی» کریجینگ عمومی» کوکریجینگ. کریجینگ با روند خارجی و رگرسیون کریجینگ 
ارزیابی شدند. برای تحلیل واریوگرافی در روش‌های کریجینگ. پنج مدل نیم‌تغیبرنمای کروی 
نمایی. «da‏ خطی دارای حد آستانه و گوسی بر داده‌های 1710 برازش شده و براساس مجموع 
مربعات» خطای کمتر و ضریب تببین بالاتر نیم تغیبرنمای مناسب انتخاب شد. تعیین مناسب‌ترین 
روش درونیابی بر اساس محاسبه شاخص‌های ریشه میانگین مربعات خطاء میانگین خطای اریب و 
میانگین انحراف مطلق روش‌های مختلف انجام شد. با بررسی شاخص‌های خطا معلوم شد روش 
کوکریجینگ با مدل‌نيم تغییرنمای گوسی دارای کمترین خطا بود و به عنوان بهترین روش برای 
مکانی‌کردن داده‌های 710 ماهانه و سالانه در خوزستان معرفی شد. در مدل‌نيم تغییرنمای گوسبی 
کوکریجینگ, نسبت بخش ساختاردار به کل تغیبرپذیری در بیشتر ماه‌ها ۱ بود که بیانگر ساختار 
مکانی قوی در تغییرات هم‌زمان متغفیرهای ETO‏ و ارتفاع است. اما دقت شود روش کوکریجینگ در 
بیشتر ماه‌ها دارای خاصیت کم برآوردی است. 


کلیدواژه‌ها: تبخیروتعرق مرجع؛ درون‌یابی» کوکریجینگ. خوزستان. 


تاریخ دریافت: ۱۳۹۱/۳۸ تاریخ تصویب:۱۳۹۱/۷/۲۱ 


ig a Gla e YE‏ . . . شماراییست‌ویکم. 
۱- مقدمه 
مهمترین عامل در تعیین نیاز آبی» برآورد دقیق ETO‏ در هر منطقه است که توزیع مکانی آن در 
مدیریت منابع آب. مطالعات بیلان آبی و برآورد نیازآبی گیاه ان اهمیت بسیاری دارد. روش‌های 
متعددی برای محاسبه ETO‏ وجود دارد که در اکثر آن‌ها از اندازه‌گیری نقطه‌ای در ایستگاه‌ها برای 
برآورد ETO‏ استفاده شده است که نتایج آنها فقط درمورد مناسب مناطق محلی کاربرد دارد. با توجه 
به اینکه بیشتر مطالعات در مقیاس منطقه‌ای وسیع صورت می گیرند» تعمیم این اطلاعات از حالت 
نقطه‌ای به منطقه‌ای ضرورت می‌یابد. با توسعه روش‌های مختلف درون‌یابی امکان برآورد مکانی 
0 منطقه‌ای در سطوح گسترده فراهم شده‌است. استان خوزستان با دارا بودن هزاران هکتار اراضی 
کشاورزی و جریان پنج رودخانه بزرگ بیشترین سهم منابع آب سطحی کشور را داراست اما 
برخلاف گستردگی "E‏ منطقه همواره مشکلات تامین آب کشاورزی بویژه در زمان 
خشکسالی‌ها مطرح بوده است. براین اساس برنامه‌ریزی و مدیریت صحیح منابع آب. تعیین نیاز آبی 
و برآورد دقیق ETO‏ ضرورت می‌بابد. هدف پژوهش حاضر تعیین مناسب‌ترین روش درونیابی برای 
برآورد مکانی و پهنه‌بندی ETO‏ در استان خوزستان می‌باشد. پژوهش حاضر به دنبال پاسخ به این 
پرسش است که «دقیق‌ترین روش درون‌یابی داده‌های ETO‏ در خوزستان کدام روش است؟» تا 


براساس آن بتوان پهنه‌بندی مکانی این پارامتر را در مقیاس سالانه و ماهانه انجام داد. 


۲- پیشینة پژوهش 

انجام شده‌است که در این تحقیق به برحی از آن‌ها اشاره شده‌است. در استان‌های فارس, بوشهر 
هرمزگان و کهگیلویه بویراحمد. سه روش زمین‌آماری کریجینگ معمولی؛ کریجینگ باقیمانده و 
کوکریجینگ برای تخمین ETo‏ ماهانه و سالانه ارزیابی شده‌اند. نتایج این بررسی‌ها نشان داده‌است 
روش کوکریجینگ به استثنای ماه‌های فروردین. اردیبهشت و شهریور بهترین روش برای برآورد 
٥‏ در این مناطق بوده‌است (نوشادی و سپاس colg‏ ۱۳۸۶: ۳۴۳). برای بررسی تغیبرات منطقه‌ای 


تبخیروتعرق de‏ در استان همدان» روش کریجینگ در مقیاس ماهانه و روش وزن‌دهی‌عکس 


سال یازدهم برآورد مکانی و پهنه‌بندی تبخیروتعرق... Yo‏ 
فاصله در مقیاس سالانه پیشنهاد شده‌است(نظری‌فر و همکاران (YA?‏ در استان تهران روش 
کریجینگ مناسب‌ترین شیوه برآورد منطقه‌ای ETO‏ بوده‌است (شهابی‌فر ۱۳۸۳). براساس پژوهشی 
که در اندلس اسپانیا انجام شد روش کریجینگ eslu‏ به عنوان بهترین روش پهنه‌بندی 1710 تعیین 
شد )242 :2008 (Vanderlinden & et al.,‏ در کشور quer‏ روش کریجینگ معمولی برای درون - 
یابی ETO‏ محاسبه‌شده به روش پنمن مانتیت برگزیده شده‌است )2005 (Gong & et al.,‏ در Oly‏ 
ببرای پیش‌بینی مکانی میانگین روزانه تبخیر و تعرق مرجع به‌دست‌آمده از معادله پنمن مانتیث» روش 
باقیمانده کریجینگ دقت بالاتری نشان‌داده‌است )251 :2005 (Mardikis & et al.,‏ در نواحی 
کوهستانی آراگون اسپانی؛ روش کوکریجینگ برای برآورد تبخیر و تعرق منطقه‌ای دارای بهترین نتیجه 
(Martinez-Cob., 1996:19) clos‏ نتایج پژوهش‌های اشاره شده OUS‏ دادند روش‌های درون - 
یابی نتایج متفاوتی در مناطق مختلف دارند براین اساس هدف پژوهش حاضر بررسی روش‌های 
مختلف درونیابی برای برآورد تبخیر و تعرق مرجع منطقه‌ای در استان خوزستان است. روش‌های 
متعددی برای تخمین ETO‏ وجود دارد و مطالعات متعددی درمورد مقایسه روش‌های مختلف برآورد 
این متغیر در داخحل و خارج از ایران انجام شده‌است. اما از بین روش‌های موجود. سازمان خواربار 
جهانی ' در نشریه شماره ۵۶ روش ترکیبی فائو پنمن مانتیث را به‌عنوان تنها روش استاندارد برای 
محاسبه نیا زآبی گیاه مرجع در تمام نقاط جهان معرفی نموده است (علیزاده و کمالی» WAY‏ 8( طبق 
آیین‌نامه مصرف بهینه آب در کشاورزی برای برآورد 1710 و نیازآبی کل کشور نیز از روش فائو 
a‏ ایت اسفاده شدءاست (وزارت کشاورزی aud lio (04 VA‏ داده‌های ETO.‏ بر فاشت‌شله از 
لایسیمتر ایستگاه هواشناسی دانشکده کشاورزی مشهد با داده‌های برآوردشده به روش‌های فائوپنمن - 
مانتیث و تشت تبخیرء نشان داد نزدیکترین جواب به داده‌های لایسیمتری از روش فائوپنمن‌مانتیث 
به‌دستآمده‌است(موسوی بایگی و همکاران ۱۳۸۸: 4۷). با بررسی روش‌های مختلف مناسب‌ترین 
روش برای تخمین ETO‏ در استان خوزستان» روش فائو پنمن‌مانتیث است(کشلولی و همکاران. 


۱۳۸۵). بنابراین برمبنای مطالعات اشاره‌شده 3 پژوهش‌های دیگر» روش فائو O^‏ مانتیث برای 


1. FAO 


برآورد نقطه‌ای ETo‏ در ایستگاه‌های مطالعاتی انتخاب شده‌است و سپس روش‌های درون‌یابی 
ارزیابی خواهند شد. 
۳-روش تحفیق 

در این پژوهش برای تعیین بهترین روش برآورد مکانی تبخیروتعرق مرجع. هفت روش 
درونیابی ازجمله وزن‌دهی‌عکس فاصله » اسپلاین, گرادیان خطی سه‌بعدی کریجینگ عمومی, 
کوکریجینگ. کریجینگ با روند خارجی" و رگرسیون-کریجینگ با یکدیگر مقایسه شدند. 
درادامه توصیف مختصری از این روش‌ها ارایه شده‌است. همچنین در این تحقیق از 
نرم‌افزارهای +05 و Minitab‏ و 9.3 Arc GIS‏ استفاده شد. 
در این مطالعه از داده‌های ۴۲ ایستگاه هواشناسی در داخل و خارج استان خوزستان (۱۳سینوپتیک, ۷ 
کلیماتولوژی و ۲۲ تبخیرسنجی) در یک دوره آماری YA‏ ساله(سال‌های ۲۰۰۹- ۱۹۸۲ میلادی) برای 
تخمین ETO‏ در محل ایستگاه‌ها استفاده شد که پراکندگی استگاه‌های مطالعاتی در شکل(۱) آمده- 
است. برای محاسبه مقادیر ETO‏ در ایستگاه‌ها به روش فائوپنمن مانتیث. نرم‌افزار کراپ وات به‌کار 
رفت. ورودی این نرم‌افزار پارامترهای اقلیمی متوسط حداقل دماء متوسط حداکثر دماء رطوبت نسبی؛ 
سرعت باد و تابش است. با توجه به اینکه در برخی ایستگاه‌های مطالعاتی» خلاء آماری وجود داشت؛ 
بنابراین پیش از تخمین تبخیروتعرق داده‌های کمبود با استفاده از آمار ایستگاه‌های دارای شرایط 
یکسان و با داده کامل تکمیل و بازسازی شدند. 


۱-۳- منطقه موردمطالعه 


محدوده مطالعاتی پژوهش استان خوزستان بود که با مساحت ۶ کبلومترمربع بین ۴۷ درجه 


و ۴۱ دقیقه تا ۵۰ درجه و YA‏ دقیقه طول شرقی و YA‏ درجه و ۵۸ دقیقه تا YY‏ درجه و ۴ دقبقه عرض 


].Inverse Distance Weighted (IDW) 
2. Kriging With External Drift 
3. Cropwat 
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شمالی واقع شده است. این منطقه دارای eil‏ خشک و نیمه خشک بوده و از لحاظ تویوگرافی به دو 


بخش کوهستانی و جلگه‌ای تقسیم می‌شود. (شکل O‏ 


شکل ۱: پراکندگی ایستگاه‌های مطالعاتی در محدوده استان خوزستان 


۱۳۹۰ نگارندگان,‎ xis 


۴- مبانی نظری 

۱-۴-کریجینگ 

کریجینگ نوعی روش تخمینی: مبتنی‌بر میانگین متحرک وزن‌دار است. برای تخمین مقادیر در 
EA‏ په ار کو د که راک یت دای کرد جم که 
۷ :۱۸۱). کریجینگ یک تابع خطی است که از مجموعه مشاهدات توزیع‌شده در همسایگی 
نقطه‌ای به‌دست می‌آید که می‌خواهيم تخمین بزنیم و تابع آن به‌صورت زیر است. 

Q) 
Z -YA.Z(x) جر ورد کته‎ dito فان سیر‎ 57. EN diodes 
m) 


4S POL : Ài مکانی مشاهده شده در نقطه 1 ام‎ gus مقدار‎ : ZX) 


به نمونه Xi‏ نسبت داده شده و بیانگر اهمیت نقطه ell‏ در برآورد است. ‏ : تعداد متغیرهای 
مکانی مشاهده‌شده‌است. 
تخمین گر کربجینگ دارای دو مرحله است. در مرحله اول شناخت و مدل‌سازی ساختار فضایی 
متغیرست که می‌توان به‌وسیله آنالیز نیم تغییرنما آن را بررسی نمود. مرحله دوم تخمین متغیر 
موردنظر را شامل می‌شود که به مرحله اول وابسته است. شرط استفاده از روش‌های کریجینگ. 
ایستا بودن متغیر می‌باشد که از طریق نیم تغییر نما قابل تشخیص است. همچنین توزیع داده‌ها 
باید به توزیع نرمال نزدیک باشد. 

۲-۴- کو کریجینگ 

درصورتی که به‌اندازه‌کافی از متغیر اصلی نمونه‌برداری نشده‌باشد و نتوان برآورد آماری را 
بادقت موردنظر انجام‌داد از روش کوکریجینگ استفاده می‌شود(33: 2007 (Majani,‏ درین‌موارد 
می‌توان با در نظر گرفتن همبستگی مکانی بین متغیراصلی و یک متغیر کمکی که از آن به اندازه 
کافی نمونه‌برداری شده‌است. تخمین را اصلاح کرد و دقت آن را بالا برد. این روش درمواقعی 
به‌کار می‌رود که ضریب همبستگی بین متغیر اصلی و کمکی. قوی و معنی‌دارست و بالاتر از 


۰/۵ می‌باشد. 


۳-۴- رگرسیون-کریجینگ 

این روش زمانی به کار می‌رود که در منطقه داده‌ها روندی وجود داشته باشد و در مرحله 
اول یک یا چند متغیر کمکی مستقل به‌وفور در منطقه موجود باشد. سپس بین متغیر وابسته و 
متغیرهای کمکی. همبستگی بالایی وجود داشته‌باشد. در این روش ابتدا با استفاده از یک رابطه 
رگرسیون خطی چندگانه بین متغیراصلی و متغیرهای کمکی. طول و عرض جغرافیایی و 
ارتفاع ایستگاه‌ها, روند موجود در داده‌ها محاسبه می‌شود سپس با کم‌کردن این میزان. از مقادیر 
متغیر اصلی. ey‏ باقیمانده در هر ایستگاه به‌دست می‌آید. سپس باقيمانده ها ۳ با روش 
کریجینگ عمومی دروذیابی می‌کنیم و «Y‏ درون‌یابی شده خطاها را به لایه مقادیر برآورد شده 
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به صورت زير است(2001:150 (Boer & et al.,‏ 


f ع(می)‎ d + Bx + Bay + ۵+ Y wz” (s) (۳ 


EI o یت تدش‎ D d Mi Gables 
است. همچنین‎ Si وزن مربوط به خطای ایستگاه‎ Wi می‌باشد و‎ So جغرافیایی و ارتفاع نقطه‎ 
است.‎ Si مقدار باقیمانده در ایستگاه‎ 2 (Si) مقدار متغیر درون‌یابی‌شده نقطه 50و‎ f (so) 

۴-۴- کریجینگ با روند خارجی 

هنگامی که داده‌های واقعی محدودی از متغیر اصلی در منطقه موجود باشد. این روش با 
استفاده از رابطه خطی بین متغیراصلی وکمکی متغیر اصلی را پیش‌بینی می کند. متغیرهای کمکی 
علاوه‌براینکه باید در محل داده‌های متغیر اصلی موجود باشند. می‌بایست در تمامی گره‌های 
شبکه پیش‌بینی نیز تعریف شده‌باشند. همچنین باید از داده‌هایی به عنوان متغیرکمکی استفاده 
نمود که همبستگی خطی بالایی با متغیر موردنظر دارند(2003:3 (Hengel & et al.,‏ 

۵-۴- روش وزن دهی عکس فاصله 

در معادله این روش وزن به صورت عکس فاصله به توان (P)‏ بیان می‌گردد که در این 
روش هرچه P‏ بزرگتر باشد بدین‌معناست که برآورد ایستگاه‌های نزدیک اهمیت بیشتری دارند. 
توانی که دارای کمترین خطا(ريشه میانگین مربع خطا) باشد. به عنوان توان بهینه انتخاب می- 
شود(2797 :2011 (Ha & et al,‏ پارامتر دیگری که در این روش تاثیر گذارست. تعداد نقاط 
همسایگی در برآورد مقدار نقاط مجهول است که در این مطالعه تعداد نقاط همسایگی 


(ابستگاه‌های مطالعاتی) YQYVUOY YU ۹ Y‏ بررسی خواهند Nue‏ 
۶-۴- روش اسپلاین 


oes ۳‏ جغرافیا و توسعة ناحیه‌ای e‏ بیست و یکم 
نمی‌توان به اطلاعات مربوط به پارامترهای توزیع azolo‏ متکی بود. اسپلاین‌ها را می‌توان یک 
تابع چندجمله‌ای‌تکه‌ای نامید. یعنی توابع پیچیده شامل قطعاتی از چندجمله‌ای‌ها با درجات 
مختلف و بین هر دو نقطه از فضا هستند که درمحل اتصالات نرم شده‌اند(مهدی‌زادی AYAN‏ 
۲ 


۷-۴-روش گرادبان خی سه بعدی 

در روش گرادیان‌های خطی فرض برین است که در منطقه موردبررسی نوعی روند در 
جهت‌های متفاوت وجود دارد که این روند یک چندجمله ای از درجه 8 فرض می‌شود. در 
روش گرادیان خطی سه‌بعدی فرض برین‌است که در راستای طول و عرض و ارتفاع منطقه 
روند uam‏ وجود دارد که با برازش یک رابطه رگرسیون خطی چندگانه بین متغیر اصلی 
اندازه گیری‌شده در ایستگاه هواشناسی و متغیرهای مستقل, می‌توان در نقاطی که متغیر اصلی 


اندازه‌گیری نشده‌است مقدار آن را تخمین زد )2007:33 (Majani,‏ . 


-A-F‏ آزمون نرمال بودن و نرمال سازی داده‌ها 

پیش از کاربرد روش‌های خانواده کریجینگ و روش‌های گرادیانی؛ باید تبعیت توزیع داده‌ها از 
توزیع نرمال مشخص‌شود. بنابراین پیش از اجرای این روش‌هاء نرمال‌بودن داده‌های تبخیروتعرق مرجع 
در سطح اطمینان ۸۹۵ بااستفاده‌از آماره اندرسون‌دارلینگ به کمک نرم‌فزار Minitab‏ انجام شد. 


۹-۴- متفیرها و شاخص‌های jos‏ 

به منظور ارزیابی دقت روش‌های درون‌یابی به‌کار گرفته‌شده برای تبخیروتعرق مرجع از 
روش اعتبارسنجی حذفی استفاده‌شد. برای ul‏ بهترین روش درون‌یابی» شاخحص‌های ريشه 
میانگین مربعات clas‏ میانگین خطای اریب" و میانگین انحراف مطلق" به کارگرفته شد که 


روابط آن‌ها در زیر T ND‏ 


1.Root Mean Square Error (RMSE) 
2. Mean Bias Error (MBE) 
3. Mean Absolute Error(« MAE) 


ze i‏ س۳ 
TEN AE Z(x,)]‏ 
ze x Z(x,)‏ ۳ 


N 
آم و لا :تعداد داده‌ها‎ Í مفدار برآوردشده متغییر‎ :Z(x;) el مشاهده‌شده متغييراً‎ EST LX) 


RMSE = 


(۳) 


MAE (۵) 


۵- یافته‌های تحقیق 

در این بخش تتایج درون‌یابی تبخبروتعرق مرجع به روش‌های مختلف در منطقه مورد مطالعه ارایه‌شد و 
در نهایت بهترین روش درولیابی این پارامت بر اساس معیارهای ارزیابی خطا مشخص‌شده‌است. 

۱-۵- نتایج نرمال سازی basala‏ 

نتایج اجرای آزمون اندرسون-دارلینگ در ماه‌های مختلف نشان‌داد که داده‌های ETO‏ در بیشتر 
موارد با سطح اطمینان ۸۹۵ از توزیع نرمال پیروی می‌نمایند جدول (V)‏ فقط ماه‌های مارس» آوریل 
می. ژوئن و دسامبر از توزیع نرمال تبعیت نکردند که با تبدیل مناسب به توزیع نرمال تبدیل شدند. 


جدول ا: نتایج تبدیل تابع توزیع احتمال داده‌های تبخیروتعرق مرجع به توزیع‌نرمال 


با زد نگارندگان» s‏ ۱۳۹ 


۳۲ مجلَة جغرافیا و توسعة ناحیه‌ای شمارة پیست و یکم 
۲-۵- نتایج روش کریجینگ عمومی 


پیش از کاربرد روش‌های خانواده کریجینگ باید ساختار مکانی داده‌ها بررسی شوند. بدین منظور ۵ 
مدل نیم تغییرنما ازجمله مدل‌های کروی, نمایی خطی, خطی که حد آستانه و گوسی داشتند بر داده ها 
برازش شده و از بین آن‌ها براساس مقدار مجموع مربعات خحطای e RSS)‏ نیم‌تغیبرنمای بهینه پایین‌تر 
انتخاب‌شد و ضریب تبیین( بالاتر گزینش‌شد. مشخصات مدل‌های نیم‌تغییرنمای بهینه ماه‌های مختلف 
در جلول (۲) آمده‌است. بررسی نسبت بخش ساعتاردار به کل تغییرپذیری نشان داد تقریباً در همه موارد 
تقش مولفه ساختاردار بسیار پیشتر از مولفه پی‌ساختارست (بیش از .07( بنابراین این سنجه دلالت بر 
ایده‌آل‌بودن ساختار مکانی داده‌های0 alala‏ در محدوده موردبررسی دارد. همچنین معلوم‌شد مدل 
تئوری کروی درتمامی‌موارد به‌عنوان مدل تثوری بهینه نیم تغییرنما درنظر گرفته‌می‌شود. دامنه تأثیر داده‌ها نیز 
بین ۲۶ تا ۸۵/۶ کیلومتر است. در شکل (۲) نمونه‌ای از نیم تغییرنمای تجربی و مدل تئوری کروی برازش 
داده‌شده برآن مشاهده می‌شود که مربوط به داده‌های تبخروتعرق ماه ژانویه است. پس از بررسی 
نیم‌تغییرنماهای تجربی و انتخاب مدل تثوری کروی به‌عنوان بهترین مدل نیم تغییرنماء درون‌یابی داده‌های 
ETO‏ ایستگاه مطالعانی به روش کریجینگ عمومی و با مدل نیم‌تغییرنمای کروی و GSH le i‏ انجام 
شد. هم‌زمان با عملیات درون‌یابی در هرماه» این روش نیز با تکنیک اعتبارسنجی حذفی ارزیابی‌شد و 
شاخص‌های خطا محاسبه‌گردید. 


حدول۲: مشخصات نیم تغییرنمای داده های تبخیرو تعرق مرجع 


^ بعات C/‏ دامنه تاثیر اثر قطعه‌ای 
n (Km) (C*C0) | a? os (RSS)‏ بش 


Vra 


1 . Residual Sums of Squares 


براورد مکانی و پهنه‌بندی تبخیروتعرق.. 


Semivariance 


47500.00 95000.00 142500.00 190000.00 


Separation Distance 


شکل ۲: نیم تغییرنمای تجربی و مدل تئوری کروی داده های تبخیروتعرق ماه ژانویه 
مأخذ: نگارندگان» ۱۳۹۰ 

۳-۵- نتایج روش کوکریجینگ 

در این مطالعه تبخیروتعرق مرجع به عنوان متغیر اصلی و ارتفاع متغیر کمکی در نظر گرفته شد. با به 
کاربستن آزمون معنی‌داری ضریب همبستگی بر این دو متفر معلوم شد در تمامی ماهها همبستگی این دو 
متغییر در سطح 14۵ معنی‌دارست. با بررسی نرمال‌بودن داده‌های ارتفاع به‌عنوان داده‌های متغیر کمکی با 
استفاده از آزمون اندرسوندارلینگ در سطح اطمینان 740 معلوم شد این داده‌ها درحد مطلوبی از توزیع نرمال 
تبعیت می‌کنند. باتوجه به اینکه داده‌های متغیر کمکی نیز باید از ساختار مکانی مطلوبی برخوردار باشند 
بنابراین پنج نوع مدل نیم‌تغییرنما بر داد‌های ارتفاع برازش شدند. آنالیز ساختاری داده‌های ارتفاع نشان داد 
آن‌ها دارای دامنه تأثیر ۳۷/۷ کیلومتر هستند و نسبت بخش ساختاردار به کل تغیبریذیری! است که نشان 
می‌دهد اثر قطعه‌ای و در نتیجه ساختار مکانی قوی داده‌های ارتفاع وجود ندارد. در روش کوکریجینگ 
تغییرات هم‌زمان دو متغیر اصلی و کمکی(نیم‌تغییرنمای متقاطع) نیز باید ساختار داشته‌باشد. بنابراین 


پارامترهای نیم‌تغییرنمای متقاطع تجربی و تئوری داده‌های تبخیروتعرق مرجع ارتفاع در ماه‌های مختلف 


Yt‏ مج جغرافیا و توسعة احیه‌ای شمارهة بیست و یکم 


[C(CCO)]‏ همه ۸۹٩ Laole‏ وا بود که نشان داد ساختار مکانی قوی در تغییرات هم‌زمان متغیرهای 


اصلی و کمکی است(جدول۳. پس از تعیین مدل‌تتوری نیم‌تغییرنمای گوسی به عنوان مدل بهینه» درون‌یابی 
داده‌های ETO‏ با متغییر کمکی ارتفاع در ایستگاه‌های مطالعاتی به روش کوکریجینگ انجام شد. 


جدول۳ : مشخصات نیم تغییرنمای متقاطع متغیر اصلی و کمک 


مجموع‌مربعات 
(RSS)u-‏ 


مأخذ: نگارندگان, ۱۳۹۰ 
۴-۵- نتایج روش رگرسیون-کریجینگ 
متغیرهای کمکی این روش طول و عرض جغرافیابی و ارتفاع ایستگاه‌ها بودند و رابطه بین متغیرهای 
کمکی و اصلی به صورت یک dal;‏ رگرسیون خطی چندگانه تعریف شد و معادله آن با نرم‌افزار 


Minitab‏ در هرماه به طور جداگانه استخراج گردید و مقادیر باقیمانده در هر ایستگاه محاسبه‌شد و در 


نهایت لایه باقیمانده ایستگاه‌ها با روش کریجینگ در نرم افزار GSE‏ درون‌یابی شد. سپس مقادیر درون- 


سال‌یاژدهمبرآورد مکانی و پهنه‌بندی تبخیروتعرقه. ٣‏ 
یابی شده در هر ایستگاه به مقادیر محاسبه‌شده توسط معادله رگرسیونی اضافه گردید و با مقادیر اصلی 
مقایسه شدند. برای ارزیابی این روش تکنیک اعتبارسنجی حذفی به‌کاررفت. 
۵-۵- نتایج روش کریجینگ با روند خارجی 
در این روش براساس رابطه عطی بین متغیر کمکی(ارتفاع) و اصلی ETO)‏ متغیر اصلی پیش‌بینی 
شد و متغیر کمکی باید در تمامی گره‌های شبکه پیش بینی موجود باشد. بدین‌منظور مختصات گره‌های 
شبکه از نرمافزار GSH‏ استخراج شد و به صورت یک aY‏ نقطه‌ای وارد نرم‌افزار ArcGIS‏ گشت. سپس 
با استفاده از مدل ارتفاعی polis DEM) ga i‏ ارتفاع متناظر با گره‌های شبکه استخراج‌شد. درنهایت 
مختصات گره‌های شبکه به همراه مقادیر ارتفاع آن‌ها؛ به عنوان متفیرکمکی برای محاسبه روش 
کریجینگ با روند خارجی وارد نرم‌افزار GSH‏ شدند. 
۶-۵- نتایج روش وزنی عکس فاصله 
در این روش توان‌های تا ۵ با تکنیک اعتبارسنجی حذفی و معیار ارزیابی خطای ریشه میانگین 
مربعات خطا بررسی شدند و توان ۳ که در بیشتر ماه‌ها دارای کمترین مقدار خطا بود به‌عنوان توان بهینه 
انتخاب‌شد. برای تعیین مناسب‌ترین تعداد نقاط همسایگی در برآورد منطقه‌ای 710 درون‌یابی داده‌ها 
برپایه تعداد ایستگاه‌های مطالعاتی ٩۰۳۶‏ ۱۲۰۱۸۰۲۱۲۴۲۷ و۳۰ و توان بهینه ۳ انجام شد. نتایج نشان داد 
که در بیشتر ماه‌ها تعداد ۶ نقطه همسایگی کمترین خطا را داشتند. 
E -۷-۵‏ روش اسپلاین 
برای درون‌یابی داده‌های تبخیرو تعرق مرجع به شیوه اسپلاین» از روش TPS‏ به وسیله 
Geostatistical Analyst Toolbox‏ از نرم افزار 9.3 ArcGIS‏ استفاده‌شد. 
۸-۵-نتایج روش گرادیان خطی سه‌بعدی 
در روش گرادیان خطی سه‌بعدی از طول و عرض جغرافیایی و ارتفاع هر ایستگاه به عنوان 
متغیرهای مستقل برای برآورد مقادیر تبخبروتعرق مرجع (متغیر وابسته) استفاده شد و بین متغیرهای 
مستقل و تبخیروتعرق‌مرجع یک رابطه رگرسیون خطی چندگانه برقرار شد. صورت کلی معادلات 
استفاده شده به شرح زیر است: 


EToj = dota xX; + 02*۲ + 0 (V) 


m‏ مج جغرافیا و توسعة احیه‌ای e‏ بیست و یکم 


:۴ : مقدار مقادیر تبخیروتعرق مرجع برآوردشده در نقطه | 

i عرض جغرافیایی نقطه‎ : Yidi طول جغرافیایی نقطه‎ : Xi 

:Zi‏ ارتفاع نقطه A1, 02,03 d‏ ,00 پارامترهای معادله 

با بررسی معادلات در ماه‌های مختلف مشخص شد در تمامی ماه‌ها ضرایب متغیرهای مستقل طول 

OX) ubi‏ عرض جغرافیایی(:) و ارتفاع C)‏ منفی است. این ننایج OUS‏ دهنده آن است که 

تبخیروتعرق مرجع از غرب به شرق و از جنوب به شمال منطقه کاهش می‌یابد و این کاهش از جنوب به 

شمال می‌تواند ناشی از تأثیر عرض جغرافیایی بر میزان تابش دریافتی در سطح زمین باشد. همچنین با 

افزایش ارتفاع ایستگاه‌ها مقدار آن کاهش می‌یابد که دلیل آن کاهش دما با افزایش ارتفاع و در نتیجه کاهش 

تبخیر از سطح زمین است. اما در این میان اثر کاهشی ارتفاع و عرض جغرافیایی بر تبخیر و تعرق مرجع در 

تمامی موارد بیشترست (به دلیل بزرگ‌تر بودن ضریب ارتفاع نسبت به ضرایب طول و عرض جنغرافیایی در 

معادلات رگرسیونی). در جدول(۴) معادلات رگرسیون خحطی سه‌بعدی ماه‌های مختلف ارایه شده‌است. 
جدول۴ : معادلات گرادیان خطی سه بعدی مربوط به ماه‌های مختلف ' 


معادله گرادیان خطی سه بعدی 
ETo = 11/3 - 0/000186 Z - 0/000001 X - 000003 Y‏ 
ETo = 18/4 - 0/000382 Z - 0/000002 X - 4 Y‏ 
ETo = 22/0 - 0/000716 Z - 0/000002 X - 0/000005 Y‏ 
ETo = 28/3 - 0/00134 Z - 0/000003 X - 0/000006 Y‏ 


ETo = 36/7 - 0/00159 Z - 0/000002 X - 0/000008 Y 
ETo = 63/5 - 0/000648 Z - 0/000008 X - 0/000015 Y 


ETo = 65/4 - 0/000040 Z - 0/000010 X - 0/000015 Y 
ETo = 60/0 - 0/000045 Z - 0/000008 X - 0/000014 Y 
ETo = 44/7 - 0/000195 Z - 0/000006 X - 0/000010 Y 
ETo = 24/6 - 0/000762 Z - 0/000001 X - 0/000006 Y 
ETo - 19/9 - 0/000562 Z - 0/000002 X - 0/000005 Y 
ETo < 12/4 - 0/000235 Z - 0/000001 X - 0/000003 Y 
ETo = 12374 - 0/204 Z - 0/00141 X - 0/00287 Y 


مأخذ: نگارندگان, ۱۳۹۰ 


۱. (26: طول جغرافیایی (متر) » ۷ : عرض جغرافیایی (متر) » 2 : ارتفاع (متر)) 


سال یازدھم برآورد مکانی و پهنه‌بندی تبخیروتعرق.. ٣۷‏ 

4-۵- انتخاب بهترین روش درون‌یابی برای منطقه‌ای نمودن تبخیرو تعرق مرجع 

برای انتخاب بهترین روش درونیابی داده‌های ۳۲0 ماهانه و سالانه در منطقه مطالعاتی 
هفت روش درونیابی ارزیابی شدند. تعیین بهترین روش مکانی کردن داده های ETo‏ ماهانه و 
سالانه بر اساس محاسبه شاخص‌های میانگین مربعات خطاء میانگین خطای اریب و میانگین 
خطای مطلق انجام پذیرفت که نتایج محاسبه شاخص‌ها به‌ترتیب در جداول )40 (۶) و (V)‏ 
موجود می‌باشد. بررسی شاخص ریشه میانگین مربعات خطا (جدول۴) در ماه‌های مختلف 
نشان داد روش‌های کریجینگ عمومی, کوکریجینگ و وزن‌دهی‌عکس فاصله نسبت‌به سایر 
روش‌ها در برآورد 81١‏ خطای کمتری دارند. نتایج این روش‌ها در بیشتر موارد نزدیک به هم 
به‌دست می آید اما بررسی شاخحص میانگین خطای مطلق (جدول۶) نشان داد در بیشتر ماه‌ها 
روش کوکریجینگ خطای کمتری نسبت‌به روش‌های دیگر دارد. بررسی شاخص میانگین 
خحطای اریبی(جدول) نشان داد روش‌های کریجینگ» کوکریجینگ و گرادیانی سه‌بعدی 
نسبت‌به ghe‏ روش‌ها دارای کمترین خطای اریب است. در شکل(۳) توزیم مکانی 
تبخیر و تعرق مرجع منطقه‌ای سالانه در استان خوزستان برای نمونه ارائه‌شده است. خطای 
ellas s (osos do dala‏ عا مل فر lass‏ اس 

جدول۵ : مقادیر ريشه میانگین مربعات خطا روش‌های مختلف درونیابی داده‌های ETo‏ ' 
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مأخذ: نگارندگان, ۱۳۹۰ 


۱. #: مقدار ريشه میانگین مربعات خطای روش منتخب در ماه مورد نظر 
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جدول ۶ : مقادیر خطای مطلق میانگین روش‌های مختلف درون‌یابی داده‌های ۲۴۲0 
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جدول ۷ : مقادیر میانگین خطاي اریبی روش‌های مختلف درون‌یابی داده‌های ۴0" 
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۱. #: مقدار میانگین خحطای مطلق روش منتخب در ماه مورد نظر 


Y‏ #: مقدار میانگین خطای اریبی روش منتخب در ماه موردنظر 
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تبخسیروتعرق مرجع سالا 
b‏ (میلی متر) 
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O 15104 - 1618.5 
7 16:56- 17366 
EE 17367-16870 5 
EN 1.5701 - 4 
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۱ 22138-22672 


شکل ۳: توزیع مکانی تبخیروتعرق مرجع منطقه‌ای سالانه در استان خوزستان- 
روش کو کریجینگ با مدل نیم تغییرنمای گوسی 
NEAR‏ نگارند گان» ۱۳۹۰ 


۶- نتیجه گیری و پیشنهادات 


با ارزیابی هفت روش درونیابی مختلف مشخص شد بهترین روش مکانی کردن و 
پهنه‌بندی داده‌های110 ماهانه و سالانه در استان خوزستان براساس شاخص‌های ريشه 
میانگین‌مربعات‌خطا و میانگین خطای‌مطلق. روش کوکریجینگ با مدل نیم‌تغییرنمای‌گوسی 
است. در مدل‌نیم تغییرنمای‌گوسی کوکریجینگ» نسبت بخش ساختاردار به کل تغییرپذیری در 
بیشتر ماه‌ها ۱ به‌دست‌آمد که بیانگر ساختار مکانی قوی در تغییرات هم‌زمان متغیرهای 110 و 
ارتفاع می‌باشد اما باید دقت‌شود کوکریجینگ در بیشتر ماه‌ها دارای خاصیت کم‌برآوردی است. 
با بررسی نقشه‌های ETO‏ ماهانه و سالانه به‌دست‌آمده معلوم شد در تمامی daol‏ حداقل ETo‏ 
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در مناطق شمال و شمال غرب استان و حداکثر آن در نواحی جنوبی به‌وقوع پیوسته‌است که 
علت آن کاهش عرض جغرافیایی از شمال به جنوب بوده است. همچنین کاهش 810 از غرب 
به شرق منطقه به‌دلیل شرایط توپوگرافی خوزستان بوده‌است که از غرب به شرق ارتفاع از 
سطح دریا افزايش یافته است و مرتفع‌ترین نقاط منطقه در نیمه شرقی قرار گرفته‌اند. با توجه 
به نتایج پژوهش حاضر که مشخص شد استفاده از متغییر کمکی در روش‌های درون‌یابی سبب 
افزايش دقت برآورد می‌گردد. پيشنهاد می‌شود در مطالعات آینده از روش‌های ترکیبی درویابی 
و متغییرهای کمکی پیشتر استفاده گردد. 

سپاسگزاری: نگارندگان برخود لازم می‌دانند از سازمان آب و برق‌خوزستان و سازمان 


هواشناسی کشور به‌دلیل در اختیارقراردادن داده‌های هواشناسی تشکر نمایند. 
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